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ABSTRACT 

Narrow dust rings observed around some young stars (e.g., HR 4796A) need to be confined. We 
present a possible explanation for the formation and confinement of such rings in optically thin cir- 
cumstellar disks, without invoking shepherding planets. If an enhancement of dust grains (e.g., due 
to a catastrophic collision) occurs somewhere in the disk, photoelectric emission from the grains can 
heat the gas to temperatures well above that of the dust. The gas orbits with super (sub)-Keplerian 
speeds inward (outward) of the associated pressure maximum. This tends to concentrate the grains 
into a narrow region. The rise in dust density leads to further heating and a stronger concentration 
of grains. A narrow dust ring forms as a result of this instability. We show that this mechanism not 
only operates around early-type stars that have high UV fluxes, but also around stars with spectral 
types as late as K. This implies that this process is generic and may have occurred during the lifetime 
of each circumstellar disk. We examine the stringent upper-limit on the H2 column density in the 
HR 4796A disk and find it to be compatible with the presence of a significant amount of hydrogen 
gas in the disk. We also compute the 01 and CII infrared line fluxes expected from various debris 
disks and show that these will be easily detectable by the upcoming Herschel mission. Herschel will 
be instrumental in detecting and characterizing gas in these disks. 
Subject headings: hydrodynamics; instabilities; circumstellar matter; infrared: stars 



1. INTRODUCTION 

Observations using infrared space telescopes, such as 
the Infrared Space Observatory (ISO), the Infrared As- 
tronomical Satellite (IRAS) and recently the Spitzer 
Space Telescope, have demonstrated that many stars 
exhibit the "Vega-phenomenon, " — excess infrared ra- 
diatio n by orbiting dust grains h eated by the central 
star ijZuckerman fe Becklinl 11993(1 . These grains are 
believed to lie in optically-thin, gas-poor disks that 
are likely the remains of more massive protop lanc- 
tary disks l)Backman fe Parescdll993HStrom fc Edwards! 
1993). Such disks are termed 'debris disks' because the 
grains have short life-times and must be continuously 
regenerated by larger bodi es. Famous examples of this 
phe nomenon includ e Vega l|Aumann et al.lll984D and 3 
Pic ([Aumanrill 19851 and references therein). iMever et al.l 
( 2006) presents the most recent review of debris disk ob- 
servations. 

When spatially resolved, many of these disks display 
peculiar morphologies, including rings, holes, clumps and 
warps. Properties of some well-known debris disks and 
their host stars are summarized in Table ^ Ring-like 
features appear to be prevalent. In this article, we fo- 
cus on the phenomenology of narrow rings, with the 
two most prominent exampl es being disks around the 
two A-type sta rs H R 4796A (l.lavawardhana et al.lll998t 
IKoerrier et alllTg^lSchueider et a.U[1 99fl>a -id Fnmal- 
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haut (Holl and et al. 2003; K alas et al. 2005; M arsh et all 
|2005|) . In the former case, a dust ring is seen in scat- 
tered light by HST at a distance of 70 A U away from the 
star with an un-resolved width < 17 AU (Sch neider et al.l 
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119991) . The dust disk in Fomalhaut was resolved by HST 
to be distributed in a narrow belt (width ~ 25 AU) at a 
distance of 135 AU from the star (|Kalas et al.ll2005fl . 

Narrow dust rings also appear around stars of later 
spectral types: the dust disks of two F-type stars, HD 
139664 (F5) and HD 181327 (F5/F6) were recently re- 
solved by HST as narrow belts, the former appears as a 
ring of width ~ 26 AU centered ~ 90 AU away from the 
star l|Kalas et alJl2006[). while the latter a width 36 AU 
at ~ 86 AU away (Schneider et al. 2006). In addition, 
model fitting using Spitzer spectrophotometry data indi- 
cate that dust around the solar- type star HD 105 may lie 
in a narrow b elt (width < 4 AU) located ~ 42 AU away 
from the star (|Mever et al.ll2004^ . 

While HR 4796A, HD 181327 and HD 105 are relatively 
young (a few Myrs), the ages for Fomalhaut and HD 
139664 are estimated to be ~ 10 8 years. 

Barring issues related to observational sensitivity lim- 
its, the presence of a clear outer edge for these rings is 
puzzling. The surface brightness of debris disks is likely 
dominated by grains just above the radiative blow-out 
size (j3 — 1/2, see eq. jH])- Even if these grains are 
produced by parent bodies on circular orbits that lie 
within a ring, radiation pressure launches them to ec- 
centric orbits upon creation (with e = (3/(1 — (3)). The 
narrow ring geometry is destroyed within a dynamical 
time and one obse rves a greatly broadened disk with no 
sharp outer edge (!Aj2£er e^ii_et_alJ|2P.Ql|; |Th£b^ji]t_ei_aL| 
| 2003UStrubbe fc Chian J2006HKrivov et al]l2006l iWvatti 
120061 also see H7.1|l . This likely explains the distribution 
of grain populations in AU Mic, f3 Pic and Vega: small 
grains are observed in a broad disk, whereas larger grains 
(detected in thermal infrared in the cases of f3 Pic and 
Vega) are observed in ring-like geometry. So why do some 
systems appear ring-like in scattered-light images? 

The most popular explanation for ring confinement in- 
vokes one or more massive shepherding planets orbit- 
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TABLE 1 

Properties of some spatially resolved debris disks 



Star 


Spectral 


L* 


Age 


X a 




Lm. 


Peak 


Ring Width 


Refs. 




Type 




(MyrJ 




(,™e) 


IT \ 


emission (AU) 


/ A T T\ 

(AUJ 




HR 4796A 


AO 


21 


8 ±2 


538 


>0.25 


5 x 10" 3 


70 


< 17 


1,2,3,4 


Vega 


AO 


60 


350 




0.003 


2 x 10~ 5 


> 86 




5,6,7,8 


Fomalhaut 


A3 


13 


100-300 




0.017 


4.6 x 10~ 5 


135 


25 


6,9,12,13 


(3 Pic 


A5 


8.7 


12i 8 4 


0.17 


0.04 


3 x 10~ 3 


120 


90° 


10,12,15 


HD 139664 


F5 


3.3 


300i™ 




> 8.7 x 10~ 4 


10~ 5 


~ 90 


26 


16 


HD 181327 


F5/F6 


3.1 


~ 10 




0.05 


2 x 10" 3 


86 


36 


17,18 


HD 105 


GO 


1.3 


30 ± 10 




- 0.03 


3.9 x 10" 4 


12 


4 


19 


HD 107146 


G2 


1.1 


80-200 


3.2 x 10" 6 


0.1-0.4 


1.2 x 10~ 3 


130 


85 


5,11 


HD 53143 


Kl 


0.7 


1000 ± 200 




1.2 x 10~ 4 


> 1.6 x 10~ 5 


~ 55 - 110 


> 55 


7,16 


e Eridani 


K2 


0.33 


850 




0.002 


8 x 10" 5 


65 


50 


12,14 



Ref erenc es. — (1) 
I20M); (6) mmUtjl 
(2001); (lll KrdSTa et al 



i et all 119981): (2llStauffer eTaT] I1995 |): ( 3 1 iGreaves et al.U20001); (41 ISchneider et alJ 119991) : ( 51 IWilliams 
20001): (71 ISong et all 1200011; (8 ) ISu e t al] ||200 5H : (9) | Barrado v N avascues e^5T ^97ll: f 101 [Z uckeonan . 
I2004D: d2llD?FbicQ et all 120 03); (13)lK alas et alJBOOa); (liTlGreave TTt^dTpOOa) : d5l IHeap eT^TT^OOol) : (16) 



12001) ; (ll)IArdila et al I2004D: (121 IDi Folco et alJ p) 0l): (13llK alas et a,l. H2005D : (14) Greaves c taU 12 
IKalas et alJ 120061) : (171 ISchneider et al.1 120061) : (181 iMannings fcljarTo^ 119981) : (191 IMever et alJ 120041) . 



a The stellar FUV (11.2 - 12.42 eV) photon flux, normalized by the average interstellar val ue of 1.2 X 10 7 / cm 2 / s IHabing|1968D , and evaluated 
at 70 AU fro m the star. The h i gh res olution spectra are provided by Kamp & Hauschildt ( Hauschildt et al. 199^), and the flux calibration is 
performed by Fernandez et al. 12006). 

^Dust masses are estimated from sub-millimeter fluxes assuming optically thin radiation, except for HD 139664 and HD 53143 for which the 
dust masses were calculated from the optical depth as described in [16]. ^^^^^^^^ 

c We adopt here the radial profile inversion by Artymowicz (private communication) based on observations from lHeap et al 1 1200(1) . 



ing at distances comparable to that of the ring. As 
grains spiral inward due to Poynting-Robertson (P-R) 
drag, the gravitational tug from the planets may tem- 
porarily stall the grains at special radii (mean-motion 
resonances, MM R), giving rise to an apparent narrow 
ring (see, e.g. IQzernqv et all 120001 IWilner et all 1 2002 : 



Quillen fe Thorndikd I2002t iKuchner fe HolmaiT 12003 : 
Moran et al.ll2004j) . The outward migration of a large 



planet may also trap g rains in MMRs, wi thout the need 
for P-R drag (see, e.g. IWvattil200l ^006). This process 
is similar to the resonance capture of Kuiper Belt ob- 
jects by Neptune. Alternatively, gravitational scattering 
by massive planets may prevent grains from migrating 
across the planet's orbit , leading to a dust distri bution 
with a sharp inner edge IjMoro-Ma rtm et aLll2005j) . 

However, there are problems with the planet hypothe- 
sis, even if we accept the existence of giant planets at such 
large radii. Firstly, grain-grain collisions typicaly induce 
radial diffusion at a rate orders of magnitude above that 
from P-R drag, thus weakening the capture in MMRs. 
The sharp edges formed by planet scattering may also 
be smeared out by collisions. Secondly, the ability of 
planets to confine grai ns that are launche d to elliptic 
orbits may be limited. IMoran et al.l l|2004|) noted that 
at eccentricities larger than that of the planet's orbit, 
grains spiralling inwar ds due to P-R drag can no longer 
be trapped in MMRs. IWvattj ((2006) further illustrated 
that even grains generated by a resonant parent popu- 
lation will become radially spread out and azimuthally 
axisymmetric, so long as they are small enough to expe- 
rience significant radiation pressure ((3 ~ 1/2). Further 
studies taking these effects into account are clearly war- 
ranted. 

Interestingly, iWvattl I^OOfift argued that the process of 
grain-grain collisions can itself act as a potential confin- 
ing mechanism for the narrow ring. For most observed 
debris disks, small dust grains collide destructively on 



timescales much shorter than that of P-R drag. Conse- 
quently, dust generated at large distances from the cen- 
tral star cannot penetrate the collisional barrier to reach 
the inner region of the system. This produces a dust 
ring with a sharply defined inner-edge. If this dust ring 
is also optically thick in the radial direction, so that ra- 
diation pressure cannot disperse the disk (see above), it 
is possible for the ring to rem ain narrow. 

iKenvon fc Bromley! l|2002') introduced an alternative 
theory for ring formation. In a planetesimal disk, run- 
away growth occurs at successively lat er times further 
away fr om the star (for an o vervie w see iNagasawa et all 
12006ft . IKenvon fc Bromlevl l|2002ft suggested that these 
newly formed protoplanets gravitationally stir the left- 
over planetesimals, inducing violent collisions that erode 
the planetesimal disk. As such, at any given time the 
disk will appear as a broad ring with a dark inner hole 
(or gap). The observed rings herald the locations of re- 
cent runaway growth and the dark holes where growth 
has c ompleted, or shadows whe re planets have yet to 
form l)Kenvon fc Bromlevl^OO^ . The dust masses and 
luminosities they obtained are com parable to those in 
the ob served debris disks (see also IKenvon fc Bromlevl 
2004b). However, if these disks are optically thin in the 
radial direction, the ring should not have a well-defined 
outer edge since the freshly produced grains will spread 
o ut in a dynamical timescale. 

iKrauss fc Wurml l)2005t) recently considered a novel ef- 
fect called photophoresis. They showed that in optically 
thin circumstellar disks with gas densities comparable to 
that in a minimum mass solar nebula (MMSN), grains 
can be pushed outward to a critical radius where the gas 
'photophoretic' force equals the gas drag. The width of 
the resulting dust ring is, as of yet, unexplored in this 
model. Even so, this model faces the following two diffi- 
culties: (1) For the critical radius to be at 70 AU, the gas 
density needs to be twice that in the MMSN. This creates 
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a problem - there can be no photophoresis effect once the 
gas has blocked out the direct starlight. At this density, 
the optical depth caused by hydrogen Rayleigh scattering 
alone is TR ay = 0.4 if integrated from 10 AU to 70 AU, 
and 24 if integrated from 1 AU to 70 AU (Rayleigh opac- 
ity k ~ 10~ 3 cm 2 /g at 4000 A. iMaver fc Dusch]ll200^ . 
This model thus requires the gas disk to be truncated in- 
ward of ~ 10 AU. (2) The assumption of twice the MMSN 
density for the observed debris disks has been ruled out 
by the observed upper limits on the column densities of 
multiple gas species (not just hydrogen) in these disks 
(|Greaves et al.ll200tt IChen fc Kamd 20041. If we assume 
a gas density of 2 x 10 -5 MMSN instead (like that in our 
model), and that grains and gas are tightly coupled, the 
critical rad ius for photophoresis sh ould be at ~ 0.2 AU 
(eq. [13] of lKrauss fc Wurml i|2005fl ). This radius is even 
smaller if grains and gas are not tightly coupled. For 
these two reasons, we believe that photophoresis cannot 
explain rings in debris d isks. 

Along a different line, iTakeuchi fc Artvmowicd 1)20011 
hereafter referred to as TA01) suggested that rings 
arise naturally in disks passing through the transitional 
phase from gas-dominated (proto-planetary) disks to 
dust-dominated (debris) disks. Residual gas in these 
disks can be rotating with either super- or sub-Keplerian 
speeds depending on the local radial pressure gradient, 
whereas dust grains have sub-Keplerian speeds that de- 
pend on grain size. The grains therefore feel either a 
head-wind or a tail-wind, forcing them to migrate ra- 
dially until they are in co- rotation with the gas. TA01 
showed that this effect causes a local concentration of 
grains; in particular, when the gas pressure experiences 
a sharp drop-off, a narrow dust ring forms locally. In 
their work this pressure drop-off is provided by an artifi- 
cial cut-off in gas density over a radial distance of order 
the disk vertical scale height. 

In this study, we propose an extension to TAOl's work. 
We lift the undesirable assumption of an artificially- 
truncated gas disk and show that narrow rings form un- 
der more general initial conditions. We assume that the 
residual gas is smoothly distributed within the disk, but 
that the dust is locally enhanced. The enhancement can 
occur over a distance of tens of AU, e.g., as the result 
of a major collisional event in the underlying planetesi- 
mal disk We relax the conventional assumption of 

temperature equality between the gas and dust compo- 
nents in the disk and show that the gas heats up due to 
the dust enhancement. We demonstrate that this gas re- 
sponse leads to the confinement of the dust grains. The 
mechanism discussed in this paper is an instability: a 
slight confinement of the dust heats the gas further which 
le ads to a stronger co nfinement . 

Kla hr fc Linl l)2005|) (hereafter KL) have proposed a 
similar mechanism for the creation of narrow rings. Their 
analytical work assumes a simple relationship between 
the gas temperature and dust density, T gas cx p^ust w hh 
(3 > 0. Our thermal analysis validates this assumption 
for a large range of gas/dust masses and stellar spectral 
types. While KL study grains of a single size and ar- 
gue that narrow rings will arise from infinitesimally small 
perturbations of a smooth dust distribution, we study a 
continuous spectrum of grain sizes and a large perturba- 
tion (a dust belt). KL also speculate that the dust ring 



"freezes" after it has formed, despite the disappearance 
of the gas as the disk evolves. This speculation needs 
to be assessed by taking into account radiation pressure 
and frequent grain collisions (see above). 

In the following, we discuss the observations related to 
dust and gas masses in known debris disk systems (Q; 
these motivate our model for the transitional disk ■ 
In we present calculations of gas heating/cooling rates 
and temperatures in such a disk. Narrow dust rings are 
shown to arise in our fiducial model, as well as in a wide 
range of parameter space (iJSJ). We discuss the best trac- 
ers for gas in these disks (©, as well as major compli- 
cations to our model (S0). We conclude in iJHl 

2. DUST AND GAS MASSES IN OBSERVED DEBRIS DISKS 
2.1. Dust Mass 

Inferred dust masses for various systems are listed in 
Tabled 0m fiducial value corresponds to that observed 
in HR 4796A. To produce ~ 0.1M© of dust (size < 1 
mm), the parent bodies (planetesimals or proto-planets) 
are likely much more massive. For comparison, the cur- 
rent Kuiper belt has a mass of ~ O.IM©, while the pri- 
mordial Kuiper belt is believed to be many times more 
massive (perhaps up to a factor of 50). 

To produce the observed infrared radiation, the dust 
disk must have a vertical scale height of 

where Lir/L* is the disk infrared to stellar optical light 
ratio (~ 5 x 10~ 3 in the case of HR 4796A (Table[TJ) and 
tv is the radial optical depth of the dust disk. In our 
fiducial model, we take H/r = 0.05. 

The vertical scale height is damped by grain-grain col- 
lisions in the collisional timescale. In our fiducial disk, 
dust travelling on inclined orbits will collide with grains 
at the mid-plane at a relative velocity of 

^ « * = 0.05. (2) 

Vkcp r 

The mean free time for collisions (T^iiision) can be esti- 
mated as 

^collision , (3) 

n s a s v± 

where a s = its 2 is the geometrical cross section of the 
grain. Only collisions between particles of similar sizes 
are relevant for the dynamics, thus n s is approximately 
the number density of dust grains in a size bin (s,2s): 

f 2s dn dust , 
n s =J s —ds, (4) 

where dndust /ds' is described by eq. 0. Tcoiiiskm, esti- 
mated at 70 AU, is plotted as a function of grain size 
in Fig. [7| (dotted line) . The collision timescale for the 
smallest grains is ~ 20 orbits in HR 4796A. 

The dust vertical velocity can also be damped by gas 
drag. In the weak coupling limit, this is described by the 
following equation, 

z = nl p z-^±z, (5) 

where z is the vertical position of the grain and vth is 
the gas thermal speed. We take /o gla in, the bulk density 
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of grains, to be 1.25 g/ cm 3 . If this process alone is im- 
portant, the vertical scale height of the dust disk decays 
with a timescale of T S ettle ~ Pgrains/Pgas^th! this is also 
plotted in Fig. 

The frequent collisions and the rapid settling by 
gas drag should result in a razor-thin debris disk. 
The observed finite thickness therefore indicate s stirring 
^Stewart fc Wetheril]lll988t IKenvon et alJH999lK 

However, giant planets are not necessarily required for 
the stirring. For instance, a Kuiper belt with individ- 
ual bodies ~ 200 km in size, self-stirred to their re- 
spective surface escape velocities, and having overlap- 
ping Hill spheres, can maintain a dust scale height of 
H/r ~ 0.05, Dust is continuously wafted up by these 
bodies. The total number of such planetesimals required 
so that their Hill spheres overlap radially and vertically 
over a belt-width of 20 AU and scale height of H/r = 0.05 
is ~ 3000 x 150, or a total mass of ~ 2M m (see also 
IKenvon fc Bromlevl200ltlOoldreich et all2004D . The co- 
pious amount of dust observed may also be produced by 
collisions between these objects. 

The observed lifetime of such a massive dust disk is 
related to the collision time of the largest particles that 
are in collisional equilibrium, which is ~ 1 Myr for 1 mm 
particles (Fig. 0). 



2.2. Gas Mass 



Except 

- It 



for 

12000 ] IPeTt 
120001 IKamp 

metal absorption 



CO 
CII 



( Rob erge et al. 
(Roberge ct al. 



some det ections of 
mt et all I2005T) CI and 

am et aTll2003l iRoberge et al]l2Q06l) and heavy 
and emi ssion lines l|01ofsson et al.l 
I200H iBrandeker et al]l2004[) . most gas observations of 
debris disks have turned up with only upper l imits 
on t h e gas column den sity (for a review, see Uural 
120041: iMever et alJ 1200^1 . These translate to gas-to- 
dust ratios of less than a few llCreaves et alJ feOOOt 



Lecavelier des Etangs et all I20QH iChen fc Kamd 2004; 



ita: 

Hollenbach et alJl2005t [Roberge et al.ll2005l also see ©• 
We set our gas-to-dust mass ratio to be in the same range. 

Gas can be dissipated either via the viscous spread- 
ing associated w ith MHD turbu l ence, or by photo- 
evaporation (e.g. lAlexander et al.l 12006). In the for- 
mer case, the disk dissipates gradually with a viscous 
timescale ~ 10 Myr (for a = 10 -2 and H/r = 0.05 at 70 
AU). HR 4796A belongs to the 8 Myr old TW Hydrae 
association, whose namesake TW Hydra is still actively 
accreting. As such, the assumption of some amount of 
gas in HR 4796A does not seem unreasonable. 

Photo-evaporation models, on the other hand, predict 
that once the inner gas disk has dissipated, which oc- 
curs at ~ 10 Myr for a 2 M^ star with an io nizing flux 
of lO^s" 1 (Table 1 of lAlexander et al.l l2~006). the outer 
disk disappears in a very short interval. If this is the 
case, our assumption of a gas disk will be difficult to jus- 
tify. However, we caution that the success of the photo- 
evaporation models depends critically on the EUV flux, 
which is uncertain even for T Tauri stars. Also, while T 
Tauri stars (on which these models are based) are likely 
magnetically active and bright in EUV, this may be less 
true for A-type or F-type stars. If one adopts instead a 
flux of 10 40 s _1 , the photo-evaporation timescale rises by 
a factor of ~ 4 and viscous evolution may dominate the 
evolution at the epoch of interest. 



3. OUR FIDUCIAL MODEL 

For our fiducial model, we consider an optically thin 
disk surrounding an HR 4796A-like star (A0V, M* = 
2.5Mp,L» = 2LL Q , T cff = 10, 000 K). We assume that 
the disk is in the so-called 'transitional phase' with dust 
mass comparable to that of the residual gas, although the 
underlying planetesimal population that produces this 
dust can have a much greater mass. Our choices for the 
masses of the gas and the dust components are justified 
in 

The mass of the microscopic dust (grains with size < 
1mm) is taken to be M^ust = O.lAf® = M gas . Dust 
grains are initially distributed within a broad (FWHM 
of 20 AU) density enhancement centered at 70 AU. The 
mid-plane dust density goes as 



PdustW = 2.7xl0- 18 exp 



(r - 70) 2 
2ct 2 



g/cm 3 , (6) 



where the standard deviation a = 20/2.355 AU. Such a 
distribution is motivated by the possibility of a major 
collisional event in the underlying planetesimal popula- 
tion. In this article we show that the broad dust dis- 
tribution described by eq. 10 will progressively narrow 
into a sharp ring. 

Initially, dust grains are assumed to follow a size dis- 
tribution of dn/ds oc s~ 4 at every radius, where s is 
the grain ra dius. This diffe rs from the classical s~ 3 - 5 
distribution l|Dohnanvilll969(K expected from collisional 
fragmentation equilibrium, but i s compatible with the re- 
sult from zodiacal dust analysis llFixsen fc DweU 2002D a s 
well as the spectral fit for HR 4796 A l|Currie et all2003jh 
With such a distribution, grains from each logarithmic 
size bin contribute comparably to the total mass, while 
the smallest grains dominate the total cross-section. Tak- 
ing s m i n = 7.7/zm (the radiative blow-out size in the HR 
4796 A system, see below), and s max = 1mm, with a grain 
bulk density of /O gra in = 1-25 g/ cm 3 , we obtain the fol- 
lowing number distribution, 



dn 



dust 



ds 



(r, s) = 1.0 x 10~ 19 



exp 



(r - 70) 2 



2<j 2 



(7) 



The radial optical depth for this dust disk is ~ 0.14. Our 
calculations take into account the associated attenuation 
of starlight. 

The dust temperature can be approximated by the 
temperature of the smallest grain at a given radius (also 
see eq. 28 of IKamp fc Bertol di 2000), 
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\AUJ 



fim 



(8) 



These grains are slightly hotter than the correspond- 
ing blackbody temperature because they are perfect ab- 
sorbers of stellar photons but inefficient emitters of their 
own thermal radiation. 

The ratio between the radiation pressure and the grav- 
ity acting on the grain reads as, 



= 



PR 



= 0.4^— — 

V s J L Q m* 



16nGM*CS Pgrain 



(9) 



where c is the speed of light and Qpr ~ 1 is the radiation 
pressure efficiency averaged over the stellar spectrum. 
The size of the smallest grain is determined by setting 
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Fig. 1. — The mid-plane gas and dust densities in our transitional 
phase disk. While the gas density (solid line) falls off as a smooth 
power-law, the dust density (dotted line) is a Gaussian centered 
at 70 AU, with a FWHM of 20 AU. The gas and dust masses 
are Af gas = M^ ust = O.lMgj and both components have the same 
vertical scale-height of H = 0.05r. We stress that although we use 
these densities for our fiducial parameters, the proposed mechanism 
works over a large range of gas and dust masses (J3- 



[3 = 0.5. 3 For an AOV star such as HR 4796A, this yields 
Smin ~ 7.7 /mi. We assume that grains move only on 
circular orbits, which limits our study to axi-symmetric 
features; see H7.1l for a discussion on eccentricity damping 
of the grains. Under this assumption, the dust orbital 
velocity (vg) is reduced from the local Keplerian value 
(ukcp = y/GM*/r) to 



v e = Vkcp (l-/?) 1/2 . 

Our gas disk has a total mass of M gas 
satisfies a mid-plane radial density profile 

/ r \ — 2.25 

3 (r)=4.0xl0- 19 (— — g/cm 3 



(10) 
0.1M ffi and 



(11) 



where r is the radial distance from the star (measured in 
AU from now on). As will become clear, the total mass 
and density profile for the gas disk matter little for our 
theory. What is relevant is the gas density at 70 AU 
and the fact that it is smoothly distributed over a broad 
region. We adopt a scale height of H/r = 0.05 for the 
gas disk, and relate surface density to mid-plane density 
simply as X = 2pH. From now on we ignore the vertical 
structure of the disk and consider only physical quanti- 
ties measured at the disk mid-plane. We also impose an 
inner and outer cut-off to the disk at 10 and 150 AU, 
respectively. These cut-offs are far from the region of in- 
terest and do not affect the dynamics. The gas density 
is plotted as a function of radius in Figure ^ 

3 Here, we define the blow-out size (s m i n ) to be the size of the 
largest grain that escapes to infinity when released on a circular 
orbit. 



The gas temperature is determined by various heat- 
ing and cooling processes (Jljl. The stellar FUV flux 
(11.2-12.42 eV) is relevant for many gas heating pro- 
cesses. Using high reso lution stellar spectra p rovided 
by Kamp & Hauschildt, iFernandez et alJ l)2006l) obtain 
a value of x = 538 for HR 4796A, where x is the stellar 
FUV flux measured at 70 AU, normalized by the Habing 
ISM field (Fa = 1.2 x 10 7 /cm 2 /s ). This is somewhat 
lar ger than the y = 460 (11.2-13.6 eV) value computed 
bv lChen fc Kamd l)2004|) from FUSE spectra. We adopt 
the latter value here. 

The most important gas heating rate in our disk model 
is photoelectric heating by dust grains. The relevant 
photon energies for this process range from W + e$ to 
13.6eV, where W is the grain work function (W ~ 8eV 
for silicate grains, and ~ 4.4 eV for carbonaceous grains, 
iWeingartner fc Draindl2001|) and $ is the electrostatic 
potential of the grain. The value of $ depends (weakly) 
on stellar spectrum, local electron density and gas tem- 
perature. It is important not to limit the photon energy 
range to the FUV range (more below). 

Gas orbiting around a star is subject to the forces of 
gravity and gas pressure. We again assume that gas 
moves only on circular orbits. The gas circular veloc- 
ity (%is) IS 



where 



g as 



n 



^kep(l - VY' 



dP_ 

dr ' 



(12) 



(13) 



Here, P = p sas kT sas / (J, g rnH is the gas pressure, and we 
take fig = 2.34 (for a mixture of molecular hydrogen and 
helium) . 

A grain that is orbiting faster than the local gas 
(vg > v gas ) experiences a head-wind, loses angular mo- 
mentum and migrates inward, and vice versa. Grains 
tend to migrate to stable orbits where they are in co- 
rotation with the gas, rj(r) — /3(s). This leads to a size 
segregation, s — s(r). This is the dynamical basis for our 
theory. 

Let the original surface brightness of a face-on disk be 
approximated by the following expression, 



J(r) 



Airr 2 



2H 



TVS 



2 C^dust , 

— - — ds 

ds 



(14) 



where L* is the stellar luminosity and H the dust scale- 
height. We set H/r = 0.05 as for the gas disk. We have 
also taken the grain albedo to be unity and assumed that 
the disk is face-on. After radial migration, grains of size 
between s and s+ds (where s > 7.7 /tm) are concentrated 
into an annulus between r and r + dr. Ignoring collisional 
evolution during this process, the dust mass within each 
ds bin is conserved. This yields a new surface brightness 
profile: 



I(r) 



L„. 



ds N s 



9 nS ~1 

47rH dr 2t 



(15) 



where N s ds is the total number of grains between s and 
s + ds integrated over the whole disk, 



dn 



dust 



ds 



2-nr'2Hdr' . 



(16) 



6 



4. GAS TEMPERATURE 

Here, we compute gas heating and cooling rates to de- 
termine the gas temperat ure profile. Much of this analy- 
sis is based on the work of lKamp fe van ZadelhofJ l|20011 
hereafter referred as KvZOl) although we do find some 
disagreements with their analysis. Except where noted, 
we adopt their expressions to calculate various rates. 

As grains migrate to their stable radii, the gas temper- 
ature profile evolves, modifying the grains' destination. 
We do not attempt to trace this behavior by solving the 
time-dependent equations self-consistently. Instead, we 
approach the problem by studying systems with dust dis- 
tributions of differing initial FWHM. 

Results at the disk mid-plane are plotted as functions 
of radius in Fig. |2 The equilibrium gas temperature is 
obtained by demanding that the local heating and cool- 
ing rates balance. The resulting temperature profile dic- 
tates the grain migration (fJSJ). Unlike KvZOl, we ignore 
heating from grain drifting (see below). 

4.1. Heating 

KvZOl conclude that the relevant gas heating processes 
in an optically thin circumstellar disk around A-stars in- 
clude collisional de-excitation of H2, photo-dissociation 
of H2, H2 formation on grains, gas-grain collisions (when 
the grain is hotter than the gas), photoelectric emission 
from dust grains and cosmic ray heating. 

While KvZOl used a fitting formula (their eq. [11]) 
to estimate the photoelectric heating rate, we opt for 
a full calculation that follows the proced ure in iDraind 
(1978) and IWeingartner & Drainei pOOll. Our results 
differ markedly from those obtained by KvZOl (see be- 
low). 

The photoelectric charging current per dust grain (of 
radius s) is 

F u 



Jo 



2 

: its e 



dv Q a b s Y{hv) 



(e<j>+W)/h 
hv-W 



hv 



f{E,hv)dE 



(17) 



where W is the work function and <j> is the charging po- 
tential of the grain. We assume that the dust disk is 
composed of 100% carbonaceous grains (W = 4.4 eV). 
Including silicon grains (W = 8eV) does not signifi- 
cantly change the photoelectric heating rate (see fJSJ. 
We adopt an overall absorption coefficient of Q a b s = 1, 
a photoelectric yield Y{hv) « (1 — W/hv)/2 and de- 
scribe the electron energy distr ibution as f(E, hv) = 
QEjjhv-W) 2 {I- E /{hv-W)) IWeingartner fc Drainei 
12001 . F v is the stellar energy flux measured at the 
grain's location, where the central star is assumed to be a 
blackbody of the appropriate temperature (10, 000 K for 
the fiducial model). F v is also attenuated by the radial 
dust optical depth. f max is the upper frequency cut-off 
of the stellar spectrum and is set at an energy of 13.6eV. 

The charging potential is obtained by equating eq. [T7j 
to the thermal electron collection current, 



J„ = Airs es e n t 



k B T, 



1 



e<f> 



(18) 



27rm e y ksTg: 
This is valid when the grains are positively charged 
(4> > 0), which is the case for the gas density of inter- 
est. We take an electron sticking coefficient of s e — 1. 



The factor of 4irs 2 arises here because the grain re- 
ceives electrons isotropically, while it intercepts photons 
only on the side facing the star (hence the factor its 2 
in eq. EI). The relevant photons for the photoelec- 
tric effect lie in the UV/FUV range, where stellar flux 
decreases (~) exponentially with wavelength. The net 
result is that the charging potential depends weakly (al- 
most logarithmically) on the electron density (n e ), gas 
temperature (T g „. s ) an d radial distance from the star (see 
Fernandez et alJl2006(l . Moreover, each grain is charged 
to a potential that is independent of its size. An easy 
result to remember is that eip ~ a few hsT* 3> fcsTgas 
over the disk. 

Short of calculating the full ionization balance in the 
transitional disk, we assume that hydrogen and carbon 
alone contribute to the electron density (n e ). Con- 
tributions from molecules and grains are irrelevant in 
our low surface density disk (at 70 AU, surface density 
~ 10~ 5 g/ cm 2 ). The ionization fraction of hydrogen is 
determined by balancin g the primary cosmic ray ioniza- 
tion rate, ~ 10 _17 s _1 (Spi tzer fc Tomaskolll9 681. with 
the rate of recombination (with recombination coefficient 
8 x IP' 11 l^fTZZ cm 3 / s. lOsterbrock fc Ferlandl 1200^1 . 
Carbon, with a first ionization potential of 11.2eV, is 
assumed to be 50% ionized. 4 Thus, 

n e = n e [B] + n e [C] 



10 



-17 



x 10 



T l/2, 
-11 1 gas ■ 



>,[H tot ] + 0.5n[C], (19) 



where n[H to t] and n[C] are the number densities of hy- 
drogen and carbon nuclei; in our model, hydrogen and 
carbon contribute comparably to the electron density. 
The associated photoelectric heating rate is 

r po = 7rs 2 / dvQ ahs Y(hv)-^ 

J(e4>+W)/h hv 



hv-W 



(E - e<j))f{E, hv)dE 



(20) 



where each ejected electron leaves the grain with a kinetic 
energy of E — e<j). 

A simple way to estimate r pc , which also clarifies the 
parameter dependences, is to set r pe ~ J e (leV/e). In 
other words, each photo-electron carries ~ leV of en- 
ergy out of the grain towards heating the gas. All other 
factors being equal, two stars that have photon fluxes 
differing by orders of magnitude can have comparable 
photoelectric heating rates (||5J). 

When grains span multiple sizes, 



r — 

x pe — 



dn 2 

ds — TTS 

ds 

hv-W 



(e<p+W)/h 

(E - e(j>)f{E, hv)dE 



Fv 

dvQ ahs Y{hv) — 



hv 



(21) 



Comparing the above heating rate against that of 
KvZOl, we find that, for the same parameters, our rates 
are typically a factor of 20 — 1000 times greater. If we 



4 T his is found to be true for the /3-Pic disk IjPern andez et alj 
2006). This fraction rises somewhat around hotter stars (like HR 
4796 A), and when the gas density (and hence the recombination 
rate) is decreased. We have not attempted to model this in detail. 
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adopt instead a cooler stellar spectrum, e.g. that of (3 
Pic, the difference increases to a factor of 10 3 — fO 5 . 
This difference arises because KvZOl used a fitting for- 
mula that only includes the FUV flux (hv > lleV). 
This is suitable for the hotter interstellar radiation field) , 
but in circumstellar environments, the star emits copi- 
ous amounts of soft UV photons (between 4.4eV and 
lOeV) which are much more relevant for grain charg- 
ing and photoelectric hea ting. Similar to our study, 
iGorti fe Hollenbachl 1)20041 appendix C4) have also con- 
sidered stellar flux between 6eV and 13.6 eV. 

KvZOl also report that the heating rate due to grain 
drift relative to the gas greatly exceeds the photoelectric 
heating rate. We find it to be unimportant in our cal- 
culations. The enhanced photoelectric heating aside, the 
drift velocity is also lower in our calculations because the 
dust and gas arc weakl y coupled. A s i milar conclusion 
was also drawn bv IGorti "fc Hollcnbach ( 200||). 

Another important heating rate comes from gas-grain 
collisions. From eq. [19] of KvZOl, 

r gg = 4.0 x 10- 12 n[H tot ] a T y%Ts 

x(T dust - T gas ) jf ns 2 ^g±ds, (22) 

where cut ~ 0.3 is the thermal accommodation coeffi- 
cient. Gas-grain collisions can contribute to either heat- 
ing or cooling, depending on the relative values of T gas 
and T dust . 

To calculate H2-related rates (H2 photo-dissociation, 
H2 formation and collisional de-excitation), we need to 
know the ratio of atomic to molecular hydrogen in the 
disk. This is determined by assuming that the H2 photo- 
dissociation and formation rates are in equilibrium. This 
is a valid approach for a star like HR 4796A, where 
the photo-dissociation timescale (~ 10 3 yrs at 70 AU) 
is much shorter than the system lifetime. The relevant 
stellar photons for H 2 photo-dissociation are in the FUV 
range (> lleV) and their flux is characterized by the 
value of x m Table ^ We include self-shielding of the 
stellar FUV flux by H2 in the disk, as well as the shield- 
ing by dust. The former accounts for the fall-off in the H2 
photo-dissociation rate after the initial rise, seen in Fig. 
121 In chemical equilibrium, H 2 formation contributes ~ 4 
times more heating than H 2 destruction. 

The conclusion that photoelectric heating is the 
dominant heating mechani sm in optically thin disks 
has also been repor t ed by I.Tonkheid et all l|2004f l and 
iKamp fc Dullemondl l)2004a|) . Furthermore, in contrast 
to the case of proto-planetary disks, the low gas density 
here implies that collisions between gas and grains can- 
not cool the gas efficiently. Instead, the gas is heated to 
temperatures well above the local dust temperature. 

4.2. Cooling 

The low gas density environment considered here is 
analogous to the surfaces of protoplanetary disks. Cool- 
ing is contributed by atomic fine structure lines from 01, 
CII, Fell and Sill, and possibly molecular transitions 
from H2, CO, H2O and OH, as well as gas-grain colli- 
sions JK^mr^fc^ ^n^a^efliofl^2fJ011|G j grti^^Iollenbach 
200l IKamp fc Dullemondl l2004at I.Tonkheid et alJ 120041: 
Kamp fc D ullemond 2004b|). Among these, we in- 
clude only the two most important processes: CII 
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10-15 " _ " 
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r (AU) 

Fig. 2. — Gas heating and cooling rates are plotted as functions 
of radius for the transitional disk around an HR 4796A-like star 
(AO). The gas is assumed to be at local thermal equilibrium. All 
cooling rates have been increased by a factor of 100 to distinguish 
them from the heating rates. At the peak of the dust enhancement 
(70 AU), photoelectric heating and OI cooling are the dominant 
heating and cooling processes, respectively. Gas-grain collisions act 
as a cooling mechanism within the vicinity of the dust enhancement 
and as a heating rate outside (not plotted) . Cosmic-rays dominate 
the heating outside the dust belt. H2-related rates are calculated 
assuming local chemical equilibrium between atomic and molecular 
hydrogen; in the absence of dust, starlight quickly photo-dissociates 
all molecular hydrogen. The asymmetric shapes of the H2-rclated 
rates around 70 AU are caused by H2 self-shielding of the stellar 
FUV flux. 



and OI. Fell (26 fim) & Sill (35 /im) can con- 
tribute a comparable amount to the cooling only 
when the gas temperature exceeds ~ 400 K. Given 
the harsh UV environment, CO, H2O and OH 
are rapidly photo-dissociated l|Kamp fc Bertoldil 120001 
ILecavelier des Etangs et aUl200ll) and are negligible for 
cooling. In comparison, H2 is abundant since it is 
continuously re-formed on grain surfaces. However, 
the lack of dipole moment and the widely spaced en- 
ergy levels preclude it from being an important coola nt 
lIHollenbach fc McKedll979t IGorti fc Hollenbachll200l . 

With a transition wavelength of 158 fim, fine struc- 
ture cooling by CII is relevant only for gas cooler than 
~ 100 K. At hotter temperatures, the fraction of CII oc- 
cupying the higher level reaches unity and the CII cool- 
ing rate saturates. OI fine-structure cooling (dominant 
line at 63 /im) kicks in at these temperatures. This rate 
saturates in turn at ~ 500 K. 

We assume standard ISM abundances for carbon and 
oxygen, n[C] = 1.4 x 10~ 4 n[H to t] and n[0] = 3.2 x 
10~ 4 n[Ht n t1, respective ly. Also, following the results in 
iFernandez et all l)2006f K we take carbon (ionization po- 
tential of 11.2eV) to be 50% ionized and oxygen (13.6eV) 
to be 100% neutral. We calculate occupation levels for 
CII and OI assuming local thermal equilibrium (LTE); 
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the Einstein coefficients and statistical weights are listed 
in Table 1 of KvZOl. When the gas density is below a 
critical value, LTE is inaccurate as the higher levels may 
not be sufficiently collisionally populated (for more de- 
tails, see KvZOl). In our model, the gas density at 70 AU 
(n[H to t] ~ 10 5 cm -3 ) lies just below the critical density 
for LTE. So our faulty LTE assumption overestimates 
the true cooling rates and underestimates the resulting 
gas temperature - a more detailed treatment will further 
strengthen the instability described here. 

We have also verified that these infrared lines are op- 
tically thin in the vertical direction. 5 

When the gas temperature exceeds the dust tempera- 
ture, gas-grain collisions act as a cooling mechanism for 
the gas. In our fiducial model, this cooling rate is sub- 
ordinate to line cooling. However, in denser disks, this 
cooling process likely dominates. 

5. THE INSTABILITY 

The gas temperature profile obtained by balancing re- 
alistic heating and cooling rates is shown in Fig. [3J Fol- 
lowing the method outlined in |]3J we compute the gas 
and dust orbital velocities, determine the locations of the 
stable orbits for different grain sizes, and calculate the 
new surface brightness profile (Fig. |3J). The latter is 
more narrowly peaked around 70 AU, with a FWHM of 
10 AU, compared to the 20 AU FWHM of the initial dis- 
tribution. More significantly, the new dust belt has both 
sharp inner and outer edges. 

We also perform the same analysis with a narrower 
initial dust distribution (FWHM = 10 AU); results are 
shown in Fig. 0] The higher initial dust density in the en- 
hancement raises the photoelectric heating rate and con- 
sequently the gas temperature. The final surface bright- 
ness is strongly peaked around 70 AU with a FWHM 
of 5 AU. Although we have not self-consistently traced 
the gas temperature evolution as the dust migrates, the 
result of this exercise suggests that the segregation and 
migration of grains acts as a positive feedback , leading to 
further evolution (also see lKlahr fe Lml l2*005). As such, 
the mechanism we discuss here is an instability. 

The profile of specific angular momentum (dashed line 
in the middle panel of Fig. |3J) decreases outwardly near 
70 AU. In this region, the gas disk will be Rayleigh un- 
stable and mix. As a result, the pressure profile will 
likely be somewhat reduced and the minimum width of 
the dust ring may be limited to ~ 10 AU. 

Two key parameters in this analysis are the gas and 
dust densities (n g , n^). The photoelectric heating rate 
scales roughly as n g rid, 6 while the gas line cooling rates 
scale as n g . So the gas temperature profile is relatively 
independent of gas density, but rises with dust density. 
In addition, the minimum size of confinable grains {(3 < 
f]max) decreases with higher dust density (also see above). 

We find that ring formation can occur over a very large 
range of dust mass: if we set an artificial limit of T max = 

5 For disks with much higher masses, these lines can become 
optically thick. This is accounted for by reducing the cooling flux 
by a factor r. 

6 The actual rate rises somewhat slower than this becaus e th e 
electron density increases slower than the gas density (cq. 1191 1. 
Also, when the gas density reaches MMSN values, grains become 
negatively charged. This leads to a Coulumb barrier that signifi- 
cantly reduces the electron collection current. 



Tdust ~ 60 K for ring formation, the minimum dust mass 
is ~ 0.005M e . The upper end for the dust mass is limited 
by dust opacity: the dust disk becomes optically thick for 
Miust > O.7M0. The range of gas mass over which this 
can occur is equally wide: the minimum gas mass is set 
by the constraint that gas is dynamically important, i.e., 
-M gas > Mdustj and the maximum gas mass approaches 
that of the MMSN value. At this value a few effects 
set in to diminish the instability including: the optical 
depth due to gas opacity, the negative charging potential 
for the grains and the increasingly important gas-grain 
collisional cooling. 

Ring formation is effective for an initial dust belt 
spread out as wide as 40 AU and is not significantly af- 
fected when a background population of dust (with 1% 
of the local gas mass) is superimposed. 

Other parameters, such as grain composition and spec- 
tral type, do not have a significant impact on our model. 
Since the work function for silicate grains is 8eV - much 
greater than the 4.4eV for graphite - photons capable 
of striking electrons out of silicate grains must be more 
energetic and therefore less numerous. As a result, the 
charging potential (<!>) for silicate grains is typically lower 
by a factor of ~ 2, which leads to a factor of ~ 2 reduc- 
tion in the photoelectric heating rate f il4.1|) . Similarly, al- 
though stars with later spectral types output drastically 
decreasing amounts of UV/FUV flux, the grain charg- 
ing potential, the electron density and the kinetic energy 
carried away by each electron only vary mildly between 
them. Indeed we find that bright, narrow rings may form 
around stars with a spectral type as late as K2. 

6. HOW TO FIND THE GAS? 

The gas mass is a critical unknown in our model. If 
gas is indeed present, what is the best way to look for it 
in debris disks? 

In Table |21 we calculate some gas observables. We as- 
sume the same fiducial parameters (e.g., gas and dust 
masses of O.IMq as in except for the characteristics 
of the central star (distance, luminosity and tempera- 
ture). We argue below that in some cases the H2 column 
densities may not be a good indicator of gas, while OI 
and CII fluxes (observable with Herschel) are excellent 
tracers. 

• H2 Column Densities 

H2 column densities have been traditionally used 
to place stringent upper limits on the gas mass. 
For instance, the non-detection of absorption in HR 
4796A establishes the following 3 — a upper-limits 
on the column densities of hydrogen: for H2 occu- 
pying the two lowest rotational levels, N[H 2! j=o] 
< 10 15 cm- 2 and N [H 2 ,j=i] < 3.7 x 10 15 cm~ 2 
llChen fc Kamdl2T)0l . 

The HR 4796 A disk plane is inclined by ~ 17° 
from our line-of-sight (iSchneider et aljfl9 9&). i.e., 
for H/r = 0.05, the line-of-sight passes ~ 5 verti- 
cal scale heights above the mid-plane . If the gas at 
70 AU is isothermal with T ~ 200 K (Fig. EJ) and 
is in vertical hydrostatic equilibrium, then the gas 
scale height is H gas ~ 23AU > H dust ~ 3.5 AU. 
If the gas at high altitudes has the same molecu- 
lar to atomic ratio as that in the mid-plane, the 
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Fig. 3. — Thermal response of the gas disk as in Fig. |2] The left panel shows gas and dust temperatures as functions of radius (in 
AU), center panel shows the corresponding gas velocities (eq. 1281 *1. normalized by the local Keplerian velocity. The dotted line marks the 
centroid (70 AU) of the dust belt; while the dashed curve indicates the profile of specific angular momentum (in arbitrary unit). Gas with 
outwardly decreasing angular momentum suffers from Rayleigh instability; we expect gas in these regions to mix rapidly and somewhat 
reduce the temperature gradient. The right panel depicts the surface brightness of a face-on disk: the dashed curve is the original surface 
brightness profile, while the solid one is the profile after the grains have migrated to their stable orbits. The surface brightness profile 
resulting from dust migration (/ ncw , right panel) is sharply concentrated around 70 AU with a FWHM of ~ 10 AU. The fluctuations seen 
are numerical artifacts. 
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Fig. 4. — Same as Fig. l3*lbut the dust is initially distributed within a narrower region (FWHM = 10 AU). Both the peak temperature 
and the maximum 77 reach higher values than those in Fig. E] The final surface brightness profile is narrowly confined between 70 and 75 
AU - the ring progressively narrows. 



line-of-sight column density of H2 will be N[H2] = 
1.2 x 10 19 cm~ 2 in our model. This is much greater 
than the observed limit. Does this rule out the 
presence of a significant amount of gas in the HR 
4796A disk? 

The gas thermal timescale at 70 AU (—nksT '/T ~ 
20 yrs for T — 200 K, where Y is the heating rate) 
is short compared to the dynamical time (~ 400 yrs 
at 70 AU), so gas rising to high altitudes will fall be- 
yond the influence of the dust and cool rapidly. The 
vertical isothermal assumption is therefore invalid 



and we do not expect much vertical expansion of 
the gas disk. Instead, H gas ~ if dust ~ 3.5 AU. This 
reduces the line-of-sight column density to N[H.2] 
~ 7 x 10 15 cm -2 . Furthermore, H2 molecules lo- 
cated well above the disk mid-plane are destroyed 
rapidly and the ratio of molecular to atomic species 
is reduced compared to that at the mid-plane. We 
therefore conclude that the observed upper limit on 
NfEfe] cannot yet exclude the presence of a signif- 
icant amount of gas in the disk. The upper limits 
on other gases (jChen fc Kamd l2004) are less con- 
straining than that of H2 and are compatible with 
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Fig. 5. — A graphical demonstration of the evolution of the dust ring, seen in scattered light. The surface brightness profiles are taken 
from Figs. I3l4l with the color white indicating brighter regions. The initial dust enhancement (left) has a FWHM of 20 AU. The gas within 
the enhancement heats up and concentrates grains into a distribution with a smaller FWHM (10 AU, center image) and sharp inner and 
outer edges. Subsequent evolution of this new distribution produces a bright ring with a FWHM of 5 AU. We have assumed that small 
grains are regenerated from larger grains after being initially lost from the system. 



TABLE 2 

LTE PREDICTIONS FOR VARIOUS DEBRIS DISKS 



Star/Spec. Type Distance Column Densities' 1 (cm 2 ) Line Fluxes at Earth b (erg/ cm 2 / s) 

(pc) N[H 2 ] N[H] OI(44.1(um) OI(63.2 /mn) OI(145.5/Ltm) CII(157.7 /an) 



HR 4796A(A0) 


67 


1.0 


X 


10 20 


1.9 


X 


10 21 


7.4 


X 


10" 


-20 


1.9 


X 


10" 


-14 


3.0 


X 


io- 


-16 


5.3 


X 


io- 


-16 


/3 Pic(A5) 


19 


2.3 


X 


10 20 


1.7 


X 


10 21 


7.7 


X 


iir 


-20 


2.1 


X 


10" 


-13 


3.1 


X 


10" 


-15 


6.0 


X 


io- 


-15 


HD 105(G0) 


40 


9.3 


X 


10 20 


2.7 


X 


10 20 


8.9 


X 


10" 


-21 


2.9 


X 


10" 


-14 


3.5 


X 


io- 


-16 


1.2 


X 


io- 


-15 


HD 107146(G2) 


28.5 


6.2 


X 


10 20 


9.0 


X 


10 20 


1.7 


X 


10" 


-20 


5.6 


X 


10" 


-14 


6.8 


X 


io- 


-16 


2.3 


X 


io- 


-15 



a Column densities for molecular and atomic hydrogen assuming that the disk is edge-on and that molecular and atomic hydrogen are 
in chemical equilibrium. The latter assumption may not be accurate for stars later than AO where the H2 photo-dissociation timcscalc 
is much longer than the system life-time. 

^The disk is optically thin in these lines for our assumed parameters. Moreover, we assume that all species are in LTE. If the number 
densities of colliding particles fall below the critical density for LTE, the lines fluxes will be greatly reduced. 



our model. 

Apart from HR 4796A, there is also a strin- 
gent upper-limit on NfHo ] in th e f3 Pic disk 
iLecavelier des Etangs et all (|2001|) : N[H 2 ] < 
10 18 cm~ 2 . This disk is edge-on, so there is not 
the same complications in interpreting the column 
densities as for HR 4796A. Most likely there is not 
a si gnificant amount of ga s in the f3 Pic disk (also 
see iFernandez et al]l200filh In corroboration with 
this conclusion, the disk does not exhibit a narrow 
dust ring. 

• OI 63.2 Flux 

In the transitional disks we study, the gas is heated 
by photoelectrons from the dust grains, and cooled 
mostly by radiation in the OI 63.2 /mi line. The 
estimated line fluxes at the Earth are listed in Table 

121 

All objects should be easily detectable by the PACS 
instrument in the upcoming Herschel mission. For 
instance, PACS should be able to achieve a 5 — a 
detection of the HR 4796A system in a mere 6 min- 
utes. With a diffraction limit of ~ 8", Herschel may 
even spatially resolve some of the close-by systems. 

Metallic gas (C, O, Fe, Na...) has been de- 
tected in the Pic disk ijBrandeker et all 12001: 



Rob erge et alll2u0^ . This gas may be produced 
during collisions of dust grains and co ntains little 
or no hydrogen ijFernandez et alJl2006D . We calcu- 
late the OI/CII fluxes expected in such a hydrogen- 
poor disk and find them to be comparable to those 
listed in Table El This is because the electron den- 
sity is not significantly reduced when hydrogen is 
absent from the gas ( ^4.1(1 . The photoelectric heat- 
ing rate remains largely unchanged from that of a 
hydrogen-rich disk. So even though the total gas 
mass of the disk is very low, there is good reasons 
to expect, and therefore to search for, emission in 
fine-structure cooling lines. 

• CII 157.7 /xm Flux 

Detection of the 157.7 /im CII line has been re- 
porte d for the j3 Pic disk (using ISO. iKamp et all 
2003) with a flux of 1.8 x 10" 13 erg/ cm 2 / s (4 - a 
detection). This is ~ 30 times greater than our 
computed value of 6.0 x 10 -15 erg/ cm 2 / s (Table 
0) and is surprising. 

Ionized carbon has been observed in absorption 
with a column densit y of N[CII] « 2.5 x 10 16 cm -2 
(Rob erge et al.l2006l) . Without knowing the actual 
gas/dust masses and distributions, We can obtain 
an upper limit to the CII flux by assuming that 
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every CII observed in absorption resides in the ex- 
cited state and spontaneously radiates. Assuming a 
vertical scale height of 0.05 and no self-absorption, 
we obtain a flux of 2.4 x 10~ 14 erg/ c m 2 / s — still 
a fact or of 10 below that reported bv lKamn et al.l 
(2003). More investigation is warranted to resolve 
this discrepancy. 

Looking towards the future, the HIFI instrument 
on Herschel, with a high frequency cut-off at 1, 910 
GHz (~ 157.0 /im) and a superb spectral resolution 
of 5 x 10 6 , may be capable of detecting this line in 
debris disks and provide us with detailed kinemat- 
ics of the gas in these disks. 

Within the temperature range relevant to this 
work, the critical electron number density required 
for LTE in the OI line is ~ 10 5 /cm 3 , while it is 
~ 10/ cm 3 in the CII line. If the disk density 
falls below these values, the expected line fluxes 
are greatly reduced from our estimates and the gas 
temperature rises. Similar values pertain for colli- 
sions with hydrogens. 

7. COMPLICATIONS 

In this section we discuss some complications ignored 
in our analysis and their impact on our conclusions. This 
include the roles of grain eccentricities, grain-grain colli- 
sions and gas dynamics. 

7.1. Eccentric Orbits of Small Grains 

Throughout our analysis, we have assumed that grains 
are on circular orbits. We examine this assumption here. 

If grains are produced from larger bodies that are on 
circular orbits, they would acquire an eccentricity at 
birth, 

P s min /r,o\ 

e = e birth = = - , (23) 

where numerically, e ~ (s/s m i n ) -5 / 3 for s ~ s m ; n . Grains 
are then launched to highly eccentric orbits with their 
apastron sorted by their sizes; the dust ring consequently 
spreads out in the dynamical timescale. Assuming a 
grain size distribution of alN s /ds oc s~", where N s is 
defined in eq. |16| . and that grains reside only at their 
apastron, r w r ap = 70 X (l + e)/(l-e) AU - 140/(1 -e) 
AU (for e ~ 1), we obtain the following surface bright- 
ness profile (eq. |T4*|): 



J(r) oc 



2-a 



ds de 



1 



140. 



r3 de dr 



-14/5+3a/5 J_ 



(24) 



This profile is insensitive to the value of a (see Fig. |5J. 
Without gas drag this broad surface brightness profile 
does not evolve with time and conflicts with the appear - 
ance of a narrow ring (also see lStrubbe fc Chiand [2006). 

The presence of gas changes this scenario: grain ec- 
centricities are typically damped by collisions with gas 
particles (which are assumed to be on circular orbits). A 
grain movi ng at a veloc ity Av relative to the gas feels a 
drag force ( Kwok 1975) 



V g asK 2 h+A« 2 ) 1/2 Av, 



(25) 



where i»th = 4/3-\/8fcBTg as /7r/iTO# is the local thermal 
speed. The relative velocity (as well as the eccentricity) 



can be damped in ~ T s orbits, where T s is the dimen- 
sionless stopping time (also called the Stokes number): 

4pgrain -^kcp 1 



XL — toil 



s ^kcp 



3p g rv t h y/l + Av 2 /v. : 



(26) 



Here f^kop is the local Keplerian frequency and t s = 
TOgrainAf/Fdrag is the dimensional stopping time, with 

^grain 

= 47r/3p gra i n s 3 . In our fiducial disk, the gas den- 
sity is so low that all grains are weakly coupled to the 
gas: T s ranges from 4 for the smallest grains to 600 for 
1 mm grains. Even so, the dimensional stopping time 
is still short compared to the collision time. So large 
grains can be assumed to be placed on circular orbits af- 
ter they are produced. Very small grains (near blow-out 
size, s > s m in), however, behave differently. 

The highly eccentric orbits (ebirth ~ 1) of these very 
small grains keep them away from the denser inner gas 
disk for most of their orbits. As a result, the effective 
stopping time is longer. Moreover, these grains encounter 
a tail wind near their apastron. For very small grains, 
this positive torque can exceed the negative torque that 
they may receive at periapstron. They are then migrated 
outward in addition to being gradually circularized. If 
the gas disk has an outer edge, the eccentricities of these 
grains may actually increase (TA01). So it is unreason- 
able to assume these small grains follow circular orbits 
with semi- major axes near 70 AU (the production site). 

For the temperature profile in Fig. |3 the large-small 
grain boundary is at ~ 3.1s mrn , while for that in Fig. 21 
the boundary is at ~ 2.4s m j n . In other words, as the gas 
temperature rises due to the increasing dust concentra- 
tion, smaller and smaller grains will be retained within 
the ring (an instability). 

In Fig. Hal we show the surface brightness profile of the 
dust ring with and without gas drag. Particles of vari- 
ous sizes are produced at 70 AU with the local Keplerian 
speed. Their dynamics then evolve under the combined 
forces of gravity, radiation pressure and gas drag (where 
applicable). In the presence of gas, the dust profile pro- 
gressively evolves into a sharply defined ring. Note that 
the local pressure maximum produced by dust heating 
is essential for containing grains - if the gas has zero 
pressure gradient or if the gas temperature equals the 
local dust temperature, all grains are gradually pushed 
outward. 

Lastly, circularization of bigger grains also reduces the 
birth eccentricity of the small grains - for instance, since 
the 20 /im grains are circularized by gas before they col- 
lide to produce the 10 /^m grains (see Fig. [7J), the birth 
eccentricity for the latter group is reduced from the value 
in eq. [231 1° 

ebirth ~ (2?) 

This moves the aforementioned large-small boundary to 
an even smaller value. 

These positive feedback effects may in the end lead 
to trapping of even the smallest grains - an impossible 
feat when no gas is present. A more complete analysis 
than that performed here is required, but the presence 
of non-circular orbits does not seem to change our basic 
conclusions about narrow-ring formation. 

7.2. Grain Collisions 
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Fig. 6. — Evolution of the surface brightness of a face-on dust 
disk due to gas drag. Here, dust particles of different sizes (s~ 4 
number distribution) are initially produced at 70 AU with the local 
Keplerian velocity. The surface brightness is plotted (in arbitrary 
unit) as a function of radius (in AU): (a) when no gas is present 
(or at time 0); (b) after undergoing gas drag for 10 4 yrs; (c) for 
10 yrs. Case (c) will be observed as a narrow ring with well- 
defined edges, confirming the results obtained assuming circular 
orbits (Fig. [3J. The adopted gas parameters are from our fiducial 
model (temperature as shown in Fig. [3J. The dot-dashed line 
indicates a r — 5 power law, confirming the result in eq. 1241 for the 
gas-free case. 



Fig. 7. — Various timescales in our fiducial model are plotted 
here as functions of grain size. T m i gI . at ; on is the radial migration 
timescale for grains to move by half the original FWHM of the 
dust belt. T co iii s ; on is the collision timescale between grains of 
comparable sizes in the belt. T se ttlc is the vertical settling time 
for a grain due to gas drag. The two dot-dashed curves repr esen t 
T c i rc and T expan< j, which are relevant for eccentric grains fl j7.ll : 
the former is the timescale for gas drag to damp the grain birth 
eccentricity, while the latter is the timescale for the periaps of small 
grains to be pushed outward of, say, 100 AU, by which distance 
they can no longer affect the narrow ring appearance. The dashed 
lines denote the estimated age for the HR 4796A system and the 
orbital period at 70 AU, respectively. The Poynting-Robertson 
drag timescale is comparable to the system lifetime even for 10 /im 
grains, and is much longer for larger grains. We set 77 = for 
simplicity when calculating T m i gra tion- 



We have so far ignored the role of grain regeneration by 
collisions. The importance of collisions can be best stud- 
ied by comparing different timescales. In Fig. [7| we plot: 
the orbital period at 70 AU; the timescale for grain colli- 
sion (T C omsion) in the initial 20 AU FWHM dust belt; the 
timescale for a grain to reach its stable orbit (Tmigration); 
the timescale for the grains to vertically settle due to gas 
drag (T sctt i ); the timescale for the grain's birth eccen- 
tricity to be halved by gas drag (T circ ); and the timescale 
for the periapstron of small grains to be expanded out- 
wards of 100 AU (T cxpan d). The latter two are relevant 
for eccentric grains. We adopt the fiducial model (Fig. 
|3J) when calculating these values. 

Numerically, we quantify T m ig ra ti on as the time it takes 
for a grain originally located at 70 AU to travel a distance 
of 10 AU (FWHM/2). The radial migration velocity is 
calculated as follows, in the weak coupling limit, the tan- 
gential movement of a grain is not significantly modified 
by the gas drag, vg = v^ cp (l — (3) 1 / 2 . The gas drag associ- 
ated with the gas-grain relative velocity removes (injects) 
angular momentum from (into) the dust grain if the gas 
moves slower (faster) than the dust, causing a grain of 
size s to migrate to a radius r where f3(s) = r](r). 7 The 

7 The Poynting-Robertson drag is insignificant in comparison to 
the gas drag. 



grains migrate radially with a velocity (cf. TA01 eq. [23]) 

1 — ?/ 



Vrad 



2 



1-/3 



1 



Vkc 



(28) 



Circular grains ( s > 300 /im) typically encounter mul- 
tiple collisions - and possibly grind down to smaller 
grains - before they reach their respective stable orbits. 
This allows small grains to be regenerated continuously 
even if they are initially lost from the system due to gas 
drag. Grains smaller than ~ 300 /im can migrate to their 
stable orbits before colliding destructively with grains of 
a similar size. 

When grain eccentricities are considered (eq. O), all 
grains greater than ~ 20 /im will be trapped in the ring 
and are circularized faster than their collision lifetime - 
their collisional progeny will likely obtain a birth eccen- 
tricity closer to that in eq. than that in eq. |23j . 
and are therefore easier to trap in the ring. Grains of 
size ~ 10 — 20 /im 8 are depleted from the ring faster than 
they can collisionally break-down (T cxpand -C T coUision ) 

8 When the gas temperature is raised further due to the enhanced 
dust concentration, these sizes are reduced to smaller values. 
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- or more relevantly, they are removed faster than they 
can be regenerated from larger grains. So these grains 
will not destroy the narrow ring appearance. Very small 
grains (7 — 8 /im) are removed more slowly. However, 
their very large apastron also means that they do not 
matter much for the narrow ring appearance. 

In summary, grain collisions do not seem to adversely 
affect our theory. 

7.3. Gas Dynamics 

We have ignored all dynamical responses of the gas to 
the heating and the migration of the dust. The gas may 
expand vertically due to the heating, and expand radially 
as the dust is concentrated. Moreover, such an optically 
thin disk most likely experiences a macroscopic viscosity 
due to MHD turbulence. How does gas dynamics impact 
the ring instability? 

In [JHl we argue that there is little vertical expansion of 
the gas within the heated regions of the disk. Because 
the gas is heated or cooled faster than the vertical sound 
crossing time, gas temperature at every altitude is deter- 
mined by the local heating and cooling. Gas high above 
the mid-plane is heated to a temperature lower than that 
at the mid-plane, since the photoelectric heating yield 
there is lower. 9 The overall vertical scale height is thus 
determined, not by the mid-plane temperature, but by 
that of the dust disk. This justifies our fiducial choice of 
setting H/r of the gas and dust components to be equal. 

As dust is concentrated into a narrow ring, the gas 
is expelled away from the ring by action-reaction, which 
lowers the local gas density and reduces the pressure gra- 
dient. This negative feedback only occurs because grains 
are being concentrated. It is a higher order effect that 
acts to limit the amplitude, but not change the sign, of 
the instability. 

Viscosity tends to erase the (angular) velocity gradient 
and dissipate the gas disk. However, if the ring is ~ 
10 AU in width, the viscous time across the ring is ~ 
0.5 Myr (taking a Shakura-Sunyaev a = 10~ 2 , H/r = 
0.05). As a result, the velocity gradient can be actively 
maintained by dust heating. Gas dissipation occurs over 
a few Myr timescale and may eventually terminate the 
ring instability. 

8. CONCLUSIONS 

We have studied the thermal and dynamical conse- 
quences of an optically thin dust disk embedded in a gas 

9 This conclusion can reverse when one considers NLTE cooling. 



disk of comparable mass. We showed that photoelectric 
emission from the dust grains heats up the gas and this 
modifies the grains' orbital motion: if the dust is pro- 
duced in a belt-like region (like that of our Kuiper belt), 
the heated gas will collect dust particles into a narrower 
region with sharp edges. This resulting dust distribution 
may be associated with the dust rings observed in var- 
ious debris disks, such as that around HR 4796A, and 
possibly HD 105. 

This mechanism operates around stars with a spectral 
type at late as K2 and is valid for gas and dust masses 
over ~ 2 orders of magnitude in range. So, if every cir- 
cumstellar disk goes through a transitional stage as we 
have described here, a narrow dust ring is expected to 
form. However, it is likely that the ring will disperse as 
the gas dissipates. 

We have restricted ourselves to consider only axi- 
symmetric rings; an interesting extension of this work 
will be to explore whether non-axisymmetric features 
may also derive from gas-dust interaction. We have also 
considered the impact of eccentricities induced by radia- 
tion pressure on small grains that are released from larger 
bodies on circular orbits. This process complicates our 
theory, but does not diminish the instability. 

We conclude that the presence of narrow dust rings 
in optically thin, transitional-stage disks, such as that 
observed around HR 4796A, does not necessarily imply 
the presence of "shepherding" planets. Instead, we argue 
that the detection of narrow dust rings may indicate the 
presence of gas. We have shown that this conclusion 
is compatible with the stringent upper limit on the H2 
column density in the HR 4796A disk. 

Lastly, most of the gas cooling occurs via the OI or CII 
infrared lines. We have computed the expected fluxes for 
various systems and point out that the upcoming Her- 
schel mission will be instrumental in observing and char- 
acterizing the gas. 
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